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ABSTRACT 


The stability of a self -gravitating, non-rotating, isothermal gas layer 
threaded by a one -dimensional equipartition magnetic field, immersed in a 
rigid isothermal layer of stars, is considered with respect to waves with motions 
in and perpendicular to the B e - g e plane, where B e and g e are the equilibrium 
magnetic and gravitational field vectors. When motions are perpendicular to 
the B e - g e plane, the magnetic field hinders gravitational instability, increasing 
the minimum length necessary to produce instability by the factor (1 + a ) 1/2 , 
where a is the ratio of magnetic pressure to gas pressure. When motions are 
in the B e - g e plane, no such simple analytical solution is found. However, 
the resulting system of equations reduces to a single fourth order differential 
equation in A0 (the perturbed gas potential) that defines an eigenvalue problem 

S 

for the marginally unstable state when the four appropriate boundary conditions 
are considered. 
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I. INTRODUCTION 


In Paper I of this series (Kellman 1972a) we investigated the equilibrium 
in the z direction (i.e. above the galactic plane) of a static, isothermal, plane- 
parallel layer of gas with equipartition magnetic and cosmic-ray components, 
immersed in an isothermal layer of stars. In Paper n (Kellman 1972b) we 
considered the gravitational stability of the gas layer with respect to plane and 
axially symmetric perturbations, neglecting the magnetic and cosmic-ray com- 
ponents, with the view toward explaining the large gas structures (1-2 kpc, 

10 7 M 0 ) observed to be the principal elements of the gaseous component of 
spiral arms in the Galaxy (McGee and Milton 1964). It is our purpose here to 
modify the stability analysis of Paper II by including a one -dimensional equi- 
partition magnetic field in the initial equilibrium state. We consider only the 
modes with motions (a) in and (b) perpendicular to the B g - g e plane, where 
B e and g e represent the equilibrium magnetic and gravitational field vectors. 

Paper IV (Kellman 1972c) will consider the effect of a combined magnetic field 
and cosmic-ray gas on the stability of the gas layer. 

Field (1970) has investigated the stability of an infinite uniform self -gravitating 
gas (no stars or cosmic-rays) against disturbances propagating both along (case a) 
and perpendicular (case b) to a uniform magnetic field B q . For case (b), he 
found that the magnetic field increased the minimum length necessary to produce 
gravitational instability by the factor (1 + 2 a ) 1 / 2 , where a is the ratio of mag- 
netic pressure to gas pressure. We will find a similar but somewhat smaller 
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modification: (1 + a) 1/2 . For case (a), Field discovered two modes. The first 

is the simple Alfven mode, and the waves travel with the Alfven velocity v = 

B /(Ait p ) 1/2 . The second is more subtle and divides into two cases, depending 
o o 

on the wavelength k . When A; < A. , where k^ = (tt <v 2 z -> /G P 0 ) 1/2 is the 
Jeans' length, the solution is stable and gives acoustic or sound waves, and in 
the limit A. < < \ , the waves travel at the sound speed. When k > k , the 

J J 

solution gives an unstable mode, the instability caused by self -gravitation of the 
gas. This is the well-known Jeans' instability. 

Parker (1966) has considered the stability of a combined magnetic field, 
cosmic-ray gas, and thermal gas against disturbances propagating along the 
direction of the magnetic field (case a). Equipartition of energy between the 
components was assumed. However, self-gravitation of the gas layer was ne- 
glected; this explains why lengths significantly smaller than typical Jeans' 
lengths resulted. 


II. STABILITY ANALYSIS 

a) Motions in a Plane Perpendicular to the B e - g Plane 

The relation between the vectors B e , g e , and k and the xyz coordinate sys- 

— ► — > . 

tern is shown in Figure 1. k is the wave propagation vector; B e and g e have 

been defined above. The case B e 1 k (case b) for an infinite uniform self- 

gravitating gas is particularly interesting (Field 1970). The longitudinal mode 

with v x = 0, v y finite propagates with phase velocity v p = <v 2 z > 1/2 (1 + 2 a) 1/2 , 
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and is called, appropriately, the magnetosonie mode. The magnetic restoring 
forces are not due to curvature of the field lines as is true of the Alfven mode, 
but to pressure gradients. The criterion for gravitational instability is X > (1 + 
2a) 1/2 \ ; the presence of a magnetic field renders gravitational instability 
more difficult. 

With these introductory remarks in mind, it seems appropriate to redefine 
the problem at hand. We seek to determine how a simple one -dimensional equi- 
partition magnetic field (along x) effects the stability criterion of waves propa- 
gating perpendicular to B e (along y) in an isothermal self -gravitating gas layer 
with infinite conductivity immersed in a rigid isothermal star layer. 

The basic equations are the continuity, momentum, hydromagnetic, Poisson, 
and heat equation, and they are written as follows: 

+/> g V-v = 0 (1) 


+ P B V - -Lb-Vb + _L VB 2 + p V0= 0 

8 d t 4 7 T 877 ^ 


( 2 ) 


— — B - Vx (v x B) = 0 
o t 


(3) 


4 77 G (P g + P„) - V 2 (0 g + 4>.) = 0 


(4) 


P 


g 


< V ?z > 


(5) 
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p, p, 0, and B are, respectively, the density, pressure, gravitational potential, 
and magnetic field strength. (v 2 tz ) is the mean square z turbulent gas velocity 
dispersion, v is the gas velocity, a first order quantity. The subscripts g and 
* denote gas and star, respectively. The perturbations in p g , p g , 0 , and B may 
be written as follows : 


^ = ^eg + A P g 

( 6 ) 

P g =P eg + A Pg = < V fz> P % 

( 7 ) 

4> - <b + A 0 

*8 ^eg T ^g 

(8) 

0 = 0 e + A 0 g - 0 eg + 0 et + A 0 g 

( 9 ) 

B = B e + A B, 

( 10 ) 


where the subscript 'e' denotes the equilibrium quantities and A denotes the 
perturbed quantities. Retaining terms only to first order in the perturbed 
quantities, equations (l)-(5) become 


3 

3 t 


A A, + 


eg 


+ /? eg V-V = 0 


( 11 ) 


<v 2 z ) VA ft 


+ A p V 0 + p + 

8 e r eg 2 t 


_V(B e -AB) +p eg VA0 g = 0 


( 12 ) 


3 

3 t 


A B + v- V B e + B e (V-v) - (B e - V) v = 0 


(13) 
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( 14 ) 


4 7r G A p - V 2 A cf> = 0, 

g g 

where we have introduced a well-known vector identity into the hydromagnetic 
equation. 

The coefficients of the system (11)-(14) are all independent of t, x, and y, 
enabling us to Fourier analyze in these variables, in which case B /Bt - n, 

B/c !x - ik x , andB /By -* ik y . However, since k = k y e y + k z e z for the case 
under study, k x = 0, and equations (11)-(14), written in component form, 
become 


nAp+ikp v+fp — — + -iL. p ) v = 0 
'g y y v e 8 3 z d z e s/ z 


(15) 


B 

1 k y <v ?*> A ^g + n ^eg v y + 1 k y A B x + i k y P eg A 0 g = 0 


(16) 


-* V L> ^~ A ^g +A ^g-J— + n ^eg V z + -T- ^ AB x + J- AB x 4~ B e + P ee A< ^e = 0 ^ 

OZ 8 8 dz 8 47732 477 dz eg 3z 8 


1 k y B e V y + 



+ n A B x = 0 


(18) 


- 4 TT G A p g 



a = 0- 


(19) 


Equations (16) and (17) are the y and z components of the momentum equation; 
equation (18) is the x component of the hydromagnetic equation. The y and z 
components of equation (13) are not written since they merely express the fact 


that AB y = AB Z = 0. 
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If we write that 


then it follows that 



f P, 


eg’ 



(20) 


( 21 ) 


since B 2 /p is assumed to be independent of z. Equation (20) is consistent 

e eg 

with the presence of a stellar component. We proceed by deriving an expression 
for d/dz 4> : 


Si *• - x; T-z (P «» + B 2 /S ”> = - ~ + ^ p„ 


eg 


<V ?z> d 

(1 + ^ dl Pe * = " <V ‘ Z> f (1 + a) ’ 


P, 


eg 


( 22 ) 


where 


a = B2/8 77 p eg <v* z > = v 2 /2 (v^ z ) . 


(23) 


Equations (20)-(23) are substituted into the system (15)-(19), with the result that 


n —? T- * + i k„ v, + ( + f\ v = 0 


e g 


y y \ 3 z 


(24) 


V . AB x A ^ g n 

— + V y + 2 a_ — + =0 


(25) 


P e% i k y <V^ z ) 


< V t 2 z> 
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J_ JL A p - f (i + a) _JL + 

P * g 3 Z g ^eg 


<V 2 > 
tz 


V 


z 


2 a 




A B. 

+ f a + 


* , 3 A ^g 


B e 3 Z (v 2 > 
t z 


= 0 


(26) 


ik v v * + yr + 


f\ 

- v + n 
2 / z 


A B 

X 

IT - 


= o 


(27) 




(v 2 2 ) g \3 Z 2 7 (V 2 z ) 


(28) 


To proceed further, we need expressions for l/B e B/Bz AB X and l/p eg B/Bz A,o 


^r/’-4 <4 - A B « ^7 (1/B -> = (A + l)4r < 29 > 


— 3— A P g = 5— (A Pg/Peg) ~ A (— 1 

p B z B B z g B 6 B z 


B + f\ (30) 


eg , 


B Z IP, 


eg 


If we define e , S , and ip by the equations 


6 = A P^Peg 
S = AB x /B e 

- A <£ g /<v 2 z > 

and employ equations (29) and (30), the system (24)-(28) reduces to 


(31) 

(32) 

(33) 
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(34) 


n 


e + i ky v y + 


fe +f ) 


v z =0 


e +. 


n 


z 


- fa ) 


1 k y < V L > 


e + n v - + 2 

<v 2 > 

t Z 


v +2a8+i//=0 


a (^ +t ) 


8 + — — - 0 = 0 
d z 


ik y v y*(-ri*j) v ’* nS = 0 


- 4 77 G p 


eg 


€ + 


< V L> 


K-i) 


= o. 


(35) 


(36) 


(37) 


(38) 


The usual way to proceed is to Fourier analyze in the remaining variable z 
and then set the determinant of the coefficients of the system (34)-(38) equal to 
zero so as to avoid the trivial solution for e , v y , v z , 8 , and 0. A dispersion 
relation between n and k y is then obtained, and each solution is called a mode. 
However, not all of the coefficients in this system are independent of z, render- 
ing Fourier analysis in z a useless exercise since the dispersion relation would 
contain unknown integrals over >/> and other of the variables. An alternate 
approach is to derive from the system (34)-(38) a single differential equation 
in the unknown 0 . With the appropriate boundary conditions the equation would 
implicitly contain the desired dispersion relation. % 
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, We have been unable to reduce the system (34)-(38) to a single differential 
equation, and therefore will restrict the analysis to the neutral or marginally 
unstable state by setting n = 0: 

e+2aS+0=O (39) 


r B , \ . / B s B , A 

- — -fa e + 2 a - — + f § + - — 0=0 

10 2 / \ B z / 3 2 


(40) 


4 " G ^ e+ f_^- k 2 U = o. 


(41) 


< V ?z> 


,B z 2 


Equations (34) and (37) have not been written explicitly; they merely express 
that v y = v^ =0, the expected result that all velocities vanish in the marginally 
unstable -state. To proceed, we differentiate equation (39) and subtract equation 
(40), with the result that 

S . e/2. < 42 > 

Combining equations (39) and (42) it easily follows that 


e - 


0 


1 + a 


(43) 


and the Poisson equation (41) thus becomes 


B 2 0 

B z 2 


4 77 G p 


eg 


[_< v ?z> (1 + “) 


-k. 2 


0 = 0 . 


(44) 
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Equation (13) in Paper II of this series (Kellman 1972b) defines the scale height 
H g of the equilibrium gas layer in the z direction to be 



<*?,> 

8 77 G Pe go ’ 


(45) 


where /? ego is the value of p at the plane z = 0. Equation (44) may therefore 
be written 


B 2 0 

d z 2 


p ( z ) 

- 66 _ ^ -k 2 


2 H 2 (1 + a) ^e go 


l_ g 


0 = 0 , 


(46) 


with the boundary conditions 


° 0 (z = 0) = 0 
o z 


(47) 


1 im 0 = 0- 

z|— *oo 


(48) 


Equation (47) results because 0 = A0 / <v 2 > is an even function of z; equation 

g t Z 

(48) results because A0 g is constrained to - 0 as | z | ■- oo. 


It is useful to recall that in Paper II we considered a problem similar in 
every respect to that considered here, except that B e = 0. Therefore, as an 
obvious check on the validity of equations (46)-(48) we require that in the limit 
B e - 0 they reduce to the analogous set of equations in Paper II, equations (18) — 
(20). This requirement is indeed satisfied. In fact, equations (46)-(48) differ 
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from equations (18)-(20) only by the factor (1 + a)=(l + B go /8 t 7 p eg < v 2 z >) in 
the second term of equation (46). Since a > 0, the depth of the equivalent 
potential well is reduced, and from our discussion in Paper II the eigenvalue 
k y decreases in magnitude. Therefore, the radius r t of the marginally unstable 
state in the symmetry plane (proportional to l/k y ) increases. This is entirely 
reasonable from a physical point of view; the presence of a magnetic field en- 
hances the difficulty of gravitational instability to result from a given disturbance. 

To calculate the radius of the marginally unstable state, we recall Figure 1 
from Paper II, where r x is plotted as a function of (v 2 z > 1 /2. We merely equate 
<v?z> 1/2 (! + a) 1/2 with <v 2 z > 1 / 2 and read off the value of r ; this is equivalent 
to multiplying r t by (1 +a ) 1/2 since r l is linearly related to <v? z > 1/2 • As 
before, we choose p ego = 1 H atom/cm 3 = 0.025 M 0 /pc 3 (Weaver 1970), P to = 
0.064 M 0 /pc3 (Luyten 1968), and <vl z ) l/2 = 18 km/sec (Woolley 1958). (v 2 > 1/2 
is allowed to vary between 1 and 20 km/sec. The results are displayed in Table 1 
where we compare r x calculated with B go = 0 fj.G, 3/xG, and 5/zG. The r x (OyuG) 
are taken from Table 1 of Paper II. We see that the magnetic field increases its 
effect as <v 2 z > 1/2 decreases. For a typical interstellar value fv 2 z > 1/2 =7.5 
km/sec, a 3/uG field increases r^OfiG) by 15%; a 5/i.p field increases r t (0 /jG) by 
44%. 

In summary, if r^O^G) is the length for which a non-magnetic, 
self-gravitating, isothermal gas layer (immersed in a rigid isothermal 
star layer) just becomes gravitationally unstable, then the introduction 
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of a one -dimensional equipartition magnetic field increases (O/j-G) by the factor 

(1 + a) 1/2 = (1 + vjj2 < V( Z » 1/2 , where the relation between B e , g e , and k is 
as shown in Figure 1. This result is similar to the modification (1 + 2 a) 1/2 
induced by the presence of a uniform one -dimensional magnetic field in an 
infinite uniform self-gravitating gas (no stars) when disturbances propagate 
across B g (Field 1970). 

b) Motions in the B e - g e Plane 

The relation between B e , g e , and k and the xyz coordinate system is shown 
in Figure 2. As in Section Ila, the basic equations are the continuity, momentum, 
hydromagnetic, Poisson, and heat equation. Only the momentum equation need 
be rewritten here, since we have chosen to write the magnetic force in a slightly 
different form: 


P g -A; v + V p - _L (V x B) x B + p V = 0. 

B a t K 4 77 B 


(49) 


Introducing the perturbations from equations (6)-(10) and retaining terms only 
to first order in the perturbed quantities, equations (1), (3), (4), (5), and (49), 
written in component form, become 


3 

Tt 


A p 

g 


+ P* 


3 

g 3 y 


v 

y 


+ 


3 d \ 

Peg Ti + dl Peg ) Vz “ 0 


(50) 
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4 77 G A p - V 2 A d) - 0 - (57) 

'e 

Equations (51), (52), and (53) are, respectively, the x, y, and z components of the 
momentum equation; equations (54), (55), and (56) are, respectively, the x, y, and 
z components of the hydromagnetic equation. Since we restrict the velocity vector 
to the y-z plane, v x = 0 and from equation (54) it follows that AB X =0. 
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The coefficients of the system (50)-(57) are all independent of t, x, and y, 
allowing us to Fourier analyze in these variables (B/Bt - n, B/Bx - ik x , B/By 
- ik ). However, since k=ke+ke,k =0 and equations (50)-(57) simplify 

y y y z z x 

to 


n A p g + i ky p 


v + 
eg y 


(- 


eg 


B z d 


eg I 


V = 0 


(58) 


i k (v? ) A p + n p v- _ A B B+ik/O A = 0 (59) 

y tz 'eg y 4 ^ Z dz e y eg g 


<vj > _JL A p + A p i_ (j> + t ip v + — 1 JL A B + -L. A B -iL B 

tz 3 Z g g dz e egz 4 77 9 z y 477 y d z e 


• B e ^ 

- 1 k v -f— A B + p — A0 = 0 

y 4 77 z ^eg z 


(60) 


( B -3T + d f l B ') v - +n 


A B y = 0 


(61) 


1 k y B e V z ~ n A B z “ 0 


(62) 


4 77 G A yO g 



(63) 


If we write as before d/dz o = fp . then it follows that d/dz B p = l/2 fB . 

eg eg e c 

since B p /p eg is assumed to be independent of z. These equations are consistent 
with the presence of a stellar component. We proceed by recalling equations 
(22), (29), and (30) for the quantities d/dz <£ e , l/B e B/Bz AB y , and l/p b/Bz Ap g . 
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Equations (58)-(63) now become 


A p 

n — + i k v 

Peg y y 



= 0 


(64) 


. . n , A.B, . A4> 

i k„ + — — v, - f a + l k 


y <„2x y B_ 


eg 


<V^ 


< v V 


(65) 


B f \ Ay0 g n _ / B f \ AB y 

^ f a] + v, + 2 a _ — + f ' 


/°eg 


B z I B 


A B -s A0 

-ik.,2 a_ — f + JL __f =0 


B B z 

e 


< v ?*> 


( 66 ) 


B f 


A B. 


iBz 2 


+ - v, + n — i = 0 


B.. 


(67) 


AB. 


1 k y v z - n -g— = 0 
e 


( 68 ) 


4 TT 


?.> \Bz2 y /<V? 2 > 


= 0. 


(69) 


As we found in Section Ila, the coefficients of the system (64)-(69) are not all 
independent of z, rendering Fourier analysis in z a useless exercise. We proceed 
as before by setting n = 0, which restricts the analysis to the marginally unstable 
state. If we make the substitutions 


e = A p / p 

'g 'eg 

(70) 

t=AB/B 

y ' e 

(71) 
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7= AB Z /B ( 


(72) 


'P = A 4> g /<v 2 tz ) , 


(73) 


equations (64)-(69) become 


ik y e-fa7+ik y '/'=0 



(~ + f)T-2ik a7+i-i|c0 
\d z I y o z 


e - (JL - kjU = o. 

<V t2^ 22 ' 


(74) 

(75) 


(76) 


Equations (64), (67), and (68) are not presented here; they merely express that 
v y = v 2 = 0, a result already expected for marginal stability. 

We proceed by attempting to obtain from the system (74)-(76) a single 
differential equation in ip , which together with the appropriate boundary conditions 
would implicitly contain the solutions for k y . It appears at first glance that we 
have three equations (74)-(76) in the four unknowns e , r , 7, and 1 p . r and 7 
are not, however, independent. From equations (67) and (68) it follows that 



and equations (74)-(76) reduce to 

ik y -fa7+ik y '/'=0 


(77) 


(78) 
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+ 


(79) 



f a] e - 


2 a 

i k 


3 2 
.3 z 2 


+ lf 


3 

3 z 




y + x— 'P = o 

a z 


4 w G P e% 

< v ?z> 



- k 2 0 = 0. 


(80) 


e and y are eliminated from the system (78)-(80), and the result is a single 
fourth order differential equation fori/; : 


1 

f 


a z« a \ f<v 2 > 

' t z 



f a f\ 
+ ~ + ~ 2 / 


B 2 

7^2 


0 


4 77 G ^eg \ dlP f 27TiQ Pe g 47 rk y G ^eg f a k y 

(v 2 z > ' */ Bz + ^ <v 2 z> f < V 2 z> ” 2 


f k y k<\ 

— + fr = °- 


(81) 


The appropriate boundary conditions to consider with equation (81) are the 
following: 


QJ , 
k QJ 

N 1 

/-N 

N 

II 

O 

N— ' 

II 

o 

(82) 

3 3 

0 (z - 0) - 0 

B z 3 

(83) 

1 im ip - 0 

(84) 

jzh® 


lim 3 ^ = 0. 

1 1 BZ2 
j z| — 00 

(85) 
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Equations (82) and (83) result because 0 = A 0 g / <v? z > is an even function of z. 
Equations (84) and (85) result because A0 g and d 2 z 2 A0 g are constrained to 
-Oas | z | - 00 . 

Equations (81)-(85) form an eigenvalue problem similar to that described 
by equations (18)-(20) of Paper II and by equations (46)-(48) of this paper. Only 
certain discrete values of k y will result in 0 that satisfy (a) equation (81) and 
(b) the boundary conditions imposed by equations (82)-(85). Each discrete k y 
corresponds to a mode of the marginally unstable state; since k « l/\ y , where 
\ is a length in the y direction, each mode is characterized by a length iV y . The 
numerical solution to this problem will be treated in a later paper of this series. 

HI. DISCUSSION 

In Paper II we discussed in some detail the observations of McGee and 
Milton (1964) relating to the existence of large-scale structure in the gaseous 
component of spiral arms in the Galaxy. Typical dimensions were observed to 
be 10 7 M @ and 1-2 Kpc. We also discussed Lin's (1970) recent proposal concerning 
the excitation of density waves in the galactic disk by classic Jeans' type in?- 
stabilities in the gaseous component of the Galaxy beyond the corotation distance. 
The inclusion of an equipartition magnetic component to the stability analysis 
should improve its applicability to both of these topics. 
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APPENDIX 


MOTIONS IN THE B e -g e PLANE (ALTERNATE DERIVATION): 

THE USE OF A VECTOR POTENTIAL 

When k and B g are constrained to lie along the y axis and v lies in the y-z 

plane (Figure 2), Parker (1966) and Field (1970) have found the vector potential 

A to be a useful quantity, where 

B = V x A. ( A1 ) 

Written in this way, the requirement V • b = 0 is automatically satisfied. The 
hydromagnetic equation thus becomes 

= 0, (A2) 

from which it follows that 

— A = vxB+VS. (A3) 

Bt 

S is an arbitrary scalar potential and will be set equal to zero. A is a useful 
quantity because it is a constant of the motion. / To see this, we introduce the 
following vector identity: 

vxg = v x (V x A) = - v • VA + V(v* A) -A-Vv - Ax (V x v). (A4) 

v*A vanishes since v is constrained to the y-z plane and A has only an x com- 
ponent. Similarly, both A-Vv and A x (V x v ) vanish. Equation (A3) therefore 
becomes 


JLa- ( ;*b) 
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— A + v ■ VA = 0, 
Bt 

(A5) 

or more simply 


~ A = 0, 

dt 

(A6) 

and thus it follows that A is a constant of the motion. 


In terms of A, the continuity, momentum, hydromagnetic, 
heat equation may be written as follows: 

Poisson, and 

4 p a + v • Vp + p V • v = 0 
g g g 

(A7) 

Vp g + Pg + ^ V " A (VA) = ° 

(A8) 

A - vx B = 0 
Bt 

(A9) 

47rG(p g + p.) - V 2 (0 g + 0.) = 0 

(A10) 

p = <v^ > p , 
tz 'g* 

(All) 

where A = A? and 

X 


-i-(VxB)xB = -l [(VxVxA)xB] = 

47 7 477 


d-(V2A),B = -L^A(VA). 

(A12) 
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Recalling the perturbations in p g , p g , <f> , and B from equations (6)-(10) and 
writing the additional equation 


A = A + AA , 

e 

equations (A7)-(A11) become 

B . 3 / B d \ n 

— Ap + p __ v + p — + — P I v =0 

Bt g eg By y y eg Bz dz eg i z 

< v t z > VA ^g + A A g W> e + Peg v + ^ [( v2 A e ) (VAA) 
+ (V 2 AA) (VA e )] + p eg VA0 g = 0 

JLaA-vxB = 0 
Bt e 


477’GAp g - V 2 A 0 g = 0. 


We have retained terms only to first order in the perturbed quantities. 


The nonvanishing component of v x B is just - B v e . In addition. 

e e z x 

may write as before that d/dz p = f p and d/dz B = 1/2 fB . Since 

eg 'eg e e 


B = V x A = — A e , 
e e dz e y 


it follows that 


B =1a 
e dz e 


(A13) 

(A14) 

(A15) 

(A16) 

(A17) 

we 

(A18) 

(A19) 
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and 


e 


_d^ 

dz 2 


A = _£ B =1 fB 
e dz e 2 e 


(A20) 


Because the coefficients of the system (A14)-(A17) are independent of x, y, and t, 
we Fourier analyze in these variables (B/Bx - ik x = 0, 3 /By ik y , B/Bt - n), 
and equations (A14)-(A17) written in component form become 


A° 

n — - + 

P, 


g ■ ik 


eg 




v = 0 

Z 


(A21) 


Ap 


£ + 1 


^eg ik (v 2 ) y 
& y ' t x- 


fa . . A( A, 

. v + AA + = 0 


(v 2 > 

t Z 


(A22) 


— - A /° g - f (! + a ) — 1£ +— v 

^eg BZ 8 Pe g (V 2 z > 




2 a 
B - 


"f B aa / 3 2 _ 2 \ a » B A( ^g 

2 9z \ Bz 2 y ) J Bz (v 2 > 

*- ' tz 


(A23) 


B v + nAA = 0 

e z 


(A24) 


477 GAp - (Jt k 2 ^ A 4> = 0. (A25) 

g \Bz 2 7 g 

We have used equation (22) for d/dz <£ e . Equations (A22) and (A23) are, re- 
spectively, the y and z components of the momentum equation; equation (A24) 
is the x component of the hydromagnetic equation. 
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into equations (A26)-(A29) and by restricting the analysis to the marginally unstable 
state (n = 0). Equations (A26)-(A29) reduce to 

e +ifLAA + <// = 0 (A32) 

B 
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(A33) 



(A34) 


e and AA are eliminated from equations (A32)-(A34), and again the result is a 
single fourth order linear differential equation for 0 : 


2 B 4 . ( S7rGp e{ 

_ 0 + I — 

f Bz 4 \ f<v^ t > 


4k 2 

y 

~T 


f a 


B 2 i/; 

Bz 2 


- 87rG k 2 p 
y eg 

f < V tz> 


AttG ip 2k 4 ’ , 

+ L!£- fak 2 - fk 2 +_L) 0 = o. 


< v tz> 


f 


(A35) 


This result may be compared with equation (81) derived above. The appropriate 
boundary conditions are expressed by equations (82)-(85). 
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TABLE 1 

RADIUS OF THE MARGINALLY UNSTABLE STATE AS A 
FUNCTION OF <v 2 > 1/2 AND B 


<v 2 > 1/2 

t Z 

(km/sec) 

r (B eo = O^tG) 
(kpc) 

r ! ( B eo = 3 ^ G ) 

(kpc) 

r (B = 
(kpc) 

2.5 

0.318 

0.67 

1.00 

5.0 

0.635 

0.87 

1.13 

7.5 

0.953 

1.10 

1.37 

10.0 

1.270 

1.40 

1.60 

15.0 

1.905 

2.03 

2.12 

20.0 

2.540 

2.61 

2.71 
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FIGURE CAPTIONS 


1. The relation between B e , g e , and k and the xyz coordinate system when 
motions are perpendicular to the B e - g e plane. 

2. The relation between B e , g e , and k and the xyz coordinate system when 
motions are in the B e - g e plane. 


28 





